Diquarks: a QCD sum rule perspective 
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We propose a phenomenological QCD sum rule with an explicit diquark field to investigate the 
essential ingredients inside the hadrons. By introducing the mass and the condensate {cj)^) for 
the diquark field as parameters in the model, we find that the sum rule works well for A, and A;,. 
This implies that these A baryons can be represented by a diquark and a quark configuration. We 
also find that there is a duality relation among the parameters (m^, (0^)) for which the sum rule is 
equally good. In the limit when (0^) = 0, we find = 400 MeV, which can be thought of as the 
constituent diquark mass. 

PACS numbers: 14.40. Gx,ll. 55. Hx,24.85.+p 



I. INTRODUCTION 



Quests for fundamental correlations in strongly inter- 
acting systems can be a clue to understanding the struc- 
ture of the physical systems. For instance, the Cooper 
pair, which is a correlated electron pair, is a key concept 
to understand many-electron systems and leads to super- 
conductivity after it condenses with other pairs. Also in 
nuclear physics, two-nucleon pair correlation can be an 
important degree of freedom to describe the nucleus, such 
as in the interacting boson model. In hadron physics, 
two-quark pair correlations, or diquarks, were introduced 
in the early days of quark model to explain systematics in 
baryon spectroscopy [l|, Over the years, experimental 
findings in the structure function, fragmentation function 
and weak decay processes pointed to the existence of a 
strong diquark correlation, especially in a definite spin- 
flavor configuration, which is the so-called good diquark 
channel 3. 

In two-color QCD with massless quarks, intcrquark in- 
teractions have exactly the same strength as q-q correla- 
tions due to the Pauli-Giirsey symmetry. This symmetry 
is explicitly broken in three-color QCD, but the remnant 
is seen in the good diquark channel, which has a resem- 
blance but less strong attraction to its counterpart of qq. 
The strong diquark correlation in the good channel can 
be derived from one gluon exchange potential or through 
the instanton-induced interaction, and is found in the 
lattice calculation in the Landau gauge Q and Coulomb 
gauge 0. It is also found in the configuration with an 
additional color source to make the whole configuration 
gauge invariant 0, However, in most hadron con- 
figurations the strong diquark correlation expressible as 
an elementary constituent has not been so evident. The 
reason why diquarks cannot be treated as an elementary 
constituent can be traced back to the original work _b 
Jaffe on the tetraquark picture for the scalar nonet 
There, the wave function for a is found to be dominantly 
composed of the tt-tt state in the meson- meson basis, 
which translates into a small diquark-antidiquark com- 
ponent in their respective good channels. This means 



that when an antiquark is present, the diquark would 
break and form a quark-antiquark pair, which is energet- 
ically more favorable than a diquark. This phenomenon 
is related to the reason we have quark condensate over 
diquark-antidiquark condensate in the QCD vacuum. 

Nonetheless, there are certain configurations where the 
good diquark correlation is the strongest and its role is 
important. The color superconducting phase appearing 
at high density is an example [l^l ■ There the strong 
correlation in the good diquark channel is responsible for 
the Cooper pairing. At the hadron level, the simplest 
configuration of this sort is the baryon composed of a di- 
quark {ud) and a spectator quark {h); namely A/i, where 
h can be s, c or 6 quarks. The reason that a diquark corre- 
lation remains strong can be understood from the color- 
spin interaction, which is known to be inversely propor- 
tional to the constituent mass of the quarks. Inside the 
A; the color-spin interaction between u — h quarks cancels 
that between d~h quarks when the diquark {ud) is in the 
1=0 channel. Moreover, the color-spin interactions with 
the heavy quark itself become smaller as h becomes heavy 
Hence the diquark will retain its strong correlation 
with small uh or dh correlations. Other configurations, 
where diquark correlations are expected to remain im- 
portant are hadrons with large angular momentum (l^ 
and exotic configurations with heavy quarks ^ISjilsj . 

In this paper, we will show that, inside the A;,, diquark 
correlation is indeed strong so that the diquark can be 
treated as an elementary constituent. We further show 
that one can extract well-defined parameters for the di- 
quark, from which one can also define the constituent 
diquark mass. 



II. FORMULATION 

To carry out our task, we propose a framework for 
descri bing the A/i based on the QCD sum rule tech- 
nique [l6l [l3l ■ In the usual study of baryon spectra in 
the QCD sum rule, one calculates a correlation function 
of baryonic currents composed of the quark fields using 
the operator product expansion . If the correlation of 
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a pair of quarks is strong, the pair can be regarded ef- 
fectively as elementary degrees of freedom. Thus, we in- 
troduce an explicit scalar-diquark field (pl with the color 
charge a in 3 and calculate the correlation function of the 
Ah current given by Ja = 4>ah'^ with h being the heavy 
quark field: 



u{q) = i J d''xe'i^T[ji,{x),.hm) 

= gff(g2)+n°(g2), 



(1) 



where q = q'^jfj,- 

The dynamics of the diquark field obeys the color 
SU(3) gauge theory: 



(2) 



The condensate parts are given as 
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where = q^-ml, D,, = q^-ml, F™" ee da;(x™(l- 
xY)l{lS.^Y with A2 = x{x ~ + (1 - x)ml + xm^, 
and 
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where the covariant derivative for the diquark is given 
by Df^ = — igt°'A'^^ and is the diquark mass. We 
incorporate Eq. © into the QCD Lagrangian. As we 
emphasized before, we will be using this Lagrangian only 
in configurations where the diquark correlations remain 
strong. One should note that the high energy limit of 
the correlation in Eq. ([1]) does not have a corresponding 
high-energy limit in real QCD, as the diquark field is 
treated here as an elementary field, which is not the case 
in reality. Nevertheless, one can use this framework to 
check the validity of the description of by confirming 
that the current Ja couples strongly to the physical Kh at 
low energy and that the correlation function reproduces 
the correct mass of A/j . 

The operator product expansion (OPE) makes it pos- 
sible to write down the correlation function in terms of 
vacuum expectation values of the elementary fields stem- 
ming from nonperturbative dynamics. The OPE of both 
the even and odd part of the correlation function ([1} can 
be calculated using the fixed-point gauge. We include the 
diquark condensate {(f)^ (j)) = {(j)"^) and perform the calcu- 
lation to leading order in up to terms of condensates 
having dimension 6. Note that since the diquark field has 
dimension 1 in the A^ current Ja = <j>''h, the correlation 
function ([T]) has mass dimension 1, which is smaller than 
the usual baryonic correlators. 

The perturbative part is as follows: 
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(3b) 

where (f = q^ — {mh ~ , <f = q^ + rn^ — and 
= 1 — Amhm^ /(f. It is worth noting that, although we 
call Eq. ^ the perturbative part that is obtained by the 
perturbative calculation of our effective model, since we 
introduce the diquark as an elementary field, these terms 
contain highly non-pcrturbative interactions in the sense 
of original QCD. 
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G°H = -^^{j{P?-'^<q^-f 

for the odd part. In each part, the first contribution to 
the gluon condensate comes from gluon operators em- 
anating from the diquark and the second contribution 
comes from the quark propagator. The factors of the 
gluon condensate proportional to (Ch) subtract the 
zero momentum part of the diquark (quark) propaga- 
tor contributing to the gluon condensate, which should 
be taken into account through the diquark (quark) con- 
densate. The parts to be subtracted can be obtained 
by the contributions from the diquark (quark) conden- 
sate after expressing the condensate in terms of the gluon 
condensate [13: 
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In fact, these equations are obtained as the leading terms 
of the heavy mass expansions. Therefore, for the A^,, 
sum rule, we can use Eq. ([6]) in Eq. (j4|) and in Eq. (O and 
obtain their sum rule without heavy quark condensate. 

Both the even and odd part of the correlation func- 
tion satisfies the following Borel transformed dispersion 
relation; 



(8) 



where B.T. stands for Borel transformation and M is 
the Borel mass. After the Borel transformation, the cor- 
relation function is expressed as an expansion in 
such that higher dimension operators have more powers 
of The Borel transformation of the fifth term of 

Eq. (U) and the last term of Eq. ([S]) proportional to 
give contributions with a factor of exp[— m^/Af^J/M^, 
which are of a much higher order in the 1 /M^ expansion. 
Thus we neglect these terms in our calculation. 

We assume that the At contribution to the correlation 
function in Eq. ([l} is given by 



n{q) = 



f{q + toa) 

— m\ + u 



^(<72-so)«PE + n&PE), (9) 



where / is a normalization of the A/i propagator, n^g° = 
lip -I- n* is the correlation function calculated by OPE 
and given in Eqs.Q, (j4|) and ([5]), and sq is the thresh- 
old parameter that represents the energy scale where the 
quark-hadron duality ansatz works in the QCD sum rule. 
Our theory will only comprise a limited phase space of 
the whole QCD. However, once we introduce the diquark 
field as an effective degree of freedom, the spectral den- 
sity at high energy can be calculated perturbatively as 
asymptotic freedom remains valid also in scalar QCD. 

In principle, the sum rules for the A;; mass can be 
obtained by the following: 
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or 



AP^M^B.T:.[e{s^~s)W{s)\ 
B.T.[0(so-s)n'(s)] 



(10) 



(11) 



where i = e or o. It should be noted however, that 11° 
is proportional to to/j except for the term proportional 
to {hh), and hence is dominated by a single term when 
rrih is small. Therefore, 11*^ should not be reliable for A. 
In such a case, one should use Eq. (|TT|) with i = e to 
calculate the mass. Moreover, as a natural extension of 
the A sum rule, we will use Eq. with i = e also for 
Ac and A;,. 

We use conventional values of the quark and gluon 
condensates and the quark mass, which were deter- 
mined in other QCD sum rules. For the case of /i = 
s, we use (fG^) = (0.33 GeV)'^, (ss) = 0.8 (gg) = 



2.0 

> 

O 1.5 



1.0 



0.5 



m=0.2, A=-0.5 
m=0.2, A=-0.0 
m=0.2, A=-1.0 
MA=1.115GeV 



2.0 

> 

O 1.5 



1.0 



1.0 1.5 2.0 2.5 
Borel mass (GeV") 



0.5 







m=0.3, A=+0.3 
m=0.3, A=-0.2 
m=0.3, A=-0.7 






MA=1.115GeV - 













1.0 



1.5 



2.0 2.5 



Borel mass (GeV ) 



FIG. 1: (Color online) Borel curves for A with several sets of 
paramters. 



-0.8(0.23GeV)^ and ^ 0.12 GeV For the 

higher dimensional condensates, we adopt the vacuum 
saturation hypothesis. That is {g^iJihY) = Kg^{hh)^ 
[20I [2l] |. , with 5 = K = 1. We will also show that the 
result with different values of k gives almost the same 
result as the sum rule has only a weak dependence on 
these operators. For the sum rules of Ac and A;,, we use 
771c = 1-47 GeV and mf, = 4.6 GeV (H], and the values 
of the condensates obtained by Eq. 

The diquark quantities to^ and {4>^(f>) arc not deter- 
mined yet and can be our model parameters to be deter- 
mined within our framework. If the diquark description 
of A/i is valid, the sum rule should work well and repro- 
duce the observed Ah mass with reasonable values of the 
diquark parameters. To parametrize the variation of the 
diquark condensate, we introduce a scaling factor A with 
respect to its value obtained in the heavy diquark mass 
limit given in Eq. ([7]) with = 0.2 GeV: 



A 



1 



96(0.2GeV)2 V 



(12) 



III. RESULTS 



A 



The plots in Fig. [T] show the Borel sum rules for the A 
mass obtained by using 11*^ in Eq. (|lll) for diquark masses 
of 0.2 and 0.3 GeV for the upper and lower figure respec- 
tively with several condensate values. One finds that for 
the upper curve, the A mass is well reproduced with a 
stable Borel curve for the diquark condensate value of 
A = —0.5, which corresponds to (0^) = — (0.04GcV)^. 
The fact that good Borel stability has been obtained im- 
plies that the correlation function of Ja with the diquark 
field strongly couples to the low-lying A baryon and thus 
the description of A in terms of the diquark degrees of 
freedom works well. 

The Borel window in this case is set between Af^ = 1.3 
and Af^ = 2.2 GeV^. The lower limit is determined by 
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FIG. 2: (Color online) Borel curves for A with diflerent thresh- 
old paramters (left) and with different k. (right). 



2.8 

2.6 

?2.4 
O 

i2.2 
^ 2.0 
1.8 





m=0.2, A=-0.5 




m=0.3, A=-0.2 




m=0.4, A= 0.0 





=2.284 GeV 
5 6 



Borel mass (GeV^) 



6.4 
6.2 
> 6.0 
O 5.8- 
S 5.6 
5.4 



5.2 
5.0 



m=0.2, A=-0.5 
m=0.3, A=-0.2 
m=0.4, A= 0.0 



iWA=5.624 GeV 



4 6 8 10 12 14 16 
Borel mass(GeV^) 



FIG. 3: (Golor online) Borel curves for Ac (left) and A;, (right) 
with the good diquark parameters. 



the OPE convergence where the contribution of the con- 
densate part is less than 30% of the perturbative part, 
while the upper limit is determined by the pole domi- 
nance in which the A pole contribution exceeds the con- 
tributions coming from the second term of Eq. (|9]) . Note 
that the present sum rule has a much wider Borel win- 
dow and better stability than the usual baryonic QCD 
sum rule. The threshold sq was set at 2 GcV^ as in a 



previous study 2j|. The sum rule is rather stable under 
a small variation in the threshold parameter leading to 
the same diquark parameters. For larger variation, one 
finds that the sum rule becomes less stable, as can be 
seen in the left figure of Fig. [2j 

Next we discuss diquark parameter dependence of the 
sum rule stability. When A becomes smaller than —0.5 
for = 0.2 GeV, the Borel stability is lost and the mass 
monotonically decreases toward smaller Borel mass. For 
larger A the curve breaks down at lower Borel mass. Simi- 
lar behavior is found for a larger diquark mass of 0.3 GeV, 
as seen in the right figure. On the other hand, this time, 
the stable curve is obtained when the diqurk condensate 
is smaller A = 0.2. One can further analyze the sum 
rule for different values of the diquark mass and diquark 
condensate and find a strong correlation between the two 
parameters. That is to say, to reproduce the physical A 
mass with a stable Borel curve, for larger (smaller) di- 
quark mass, we need smaller (larger) diquark condensate. 
In general, obtaining a stable curve should have enough 
information to determine a unique pair of the parame- 
ters, since it is not trivial at all that the sum rule works 
well with the Borel stability. Nevertheless, in the present 
case, we find a duality-type relation among the two pa- 
rameters. This means that the A mass may be essentially 
controlled by a single physical quantity, such as a "con- 
stituent" diquark mass. In the limit when the diquark 
condensate is zero (A = 0), one finds that the necessary 
diquark mass is around 0.4 GeV, which corresponds to 
the constituent diquark mass. 

Finally, we show the sum rule with (X,m,p) = 
(—0.5, 02GeV) but with different k. As can be seen in 
the right figure of Fig. [51 the dashed and dot-dashed lines 
are almost indistinguishable from the solid line, suggest- 



ing a very weak dependence on the higher dimensional 
operators. 



B. Ac and A^ 

We also calculate the Ac and A;, masses using the good 
diquark parameters obtained in the A sum rule, we set 
So —7 and 40 GeV^, respectively. Figure [3] shows the 
Borel curves for Ac and A;, obtained with the sum rule 
given in Eq. (fTOl) . As can be seen in the figure, all sets of 
parameters (A, m^) that worked well for A reproduce the 
observed Ac very well. Although the best Borel curve 
slightly underestimates the mass of At,, overall our re- 
sults imply that the diquark description works well again, 
suggesting that, in these A/j baryons, the diquark config- 
uration represents important degrees of freedom having 
universal properties. 

In Fig. m we show the stability in the choice of the con- 
tinuum threshold by plotting the Borel curve for different 
So values. Here, we plot the graphs for one set of the ac- 
cepted value of (A, m^) = (—0.5, 0.2GeV). As can be seen 
from the figure, we obtain a reasonable mass of Ac when 
choosing the value of sq which gives the most stable Borel 
curve. Unfortunately, the threshold parameter that gives 
the most stable Borel curve slightly underestimates the 
mass of Ab. We believe that the large difference in the 
bottom quark mass requires the use of different renormal- 
ization points for the vacuum expectation values, which 
can be a topic of future study. 



C. Summary and discussions 

The calculation of the correlation function as given in 
Eq. ([T]) means performing an OPE, assuming that the di- 
quark structure is intact even when the relative momen- 
tum between the diquark effective field and the heavy 
quark is large. This is a simplification, but since we are 
interested only in the structure of the physical A, all we 
need is that the approximation be valid for the physical 
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A, Ac and Af,, as has been shown through our QCD sum 
rule analysis. 

On the other hand, when we are analyzing states that 
have another light quark or an antiquark, the diquark de- 
scription is expected to fail. For example, in the proton, 
the ud diquark, despite its strong correlation, will tunnel 
into a new du diquark with the spectator u quark. In our 
language, we would need an effective (j) — u ~ d coupling 
to include this process. The actual application of our 
sum rule to the nucleon case obtained by taking to 
zero gives a nucleon mass above 1 GeV without a stable 
Borel curve. In that sense the minimal Lagrangian for 
the diquark given in Eq. (2) is unique for describing A/j. 

We have also analyzed the sum rule with a = 4>^4> 
current. In the scalar nonet picture, the a meson is 
composed of four quarks, but will have a strong quark- 
antiquark component. So we expect the coupling of our 
current with the physical a to be very small. Indeed, we 
find that the sum rule obtained does not reproduce the 
a mass around 400 to 600 MeV. Instead we find a mono- 
tonically decreasing Borel curve where the mass ranges 
from 1400 to 1000 GeV. The exphcit treatment of four 
quark fields in the QCD sum rule has reproduced a lower 
mass with a stable curve [1^ . 

In conclusion, we have introduced an explicit diquark 
degree of freedom and analyzed the sum rule for A, Ac 
and Af,. We find that we can reproduce the masses of 
these states well, suggesting that these states can be rep- 
resented by a strong diquark and a spectator quark. We 
further find that there is a duality between diquark con- 



densate and diquark mass from which we can estimate 
the constituent diquark mass to be around 400 MeV. We 
expect that such an explicit diquark degree of freedom 
can be an effective method to describe configurations 
where diquark correlation is expected to remain strong. 
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